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Abstract 

The spatial search problem consists in minimizing the number of steps re- 
quired to find a given site in a network, under the restriction that only oracle 
queries or translations to neighboring sites are allowed. We propose a quantum 
algorithm for the spatial search problem on a triangular lattice with TV sites 
and torus-like boundary conditions. The proposed algorithm is a special case of 
the general framework for abstract search proposed by Ambainis, Kempe and 
Rivosh |Kempe et al. 2005] (AKR) and Tulsi |Tulsi 2008| . applied to a trian- 
gular network. The AKR- Tulsi formalism was employed to show that the time 
complexity of the quantum search on the triangular lattice is 0(\/N\ogN). 

1 Introduction 

The spatial quantum search problem consists of using local unitary operations to 
search for one (or more) nodes within a set of N spatially arranged sites with an 
implicit notion of distance between them. The search nodes are identified by the 
non-zero values of a binary function (the oracle), as usual. The spatial search prob- 
lem [Ambainis and Aaronson 2003] incorporates the restriction that, in one step, one 
can cither query the oracle at the current site or advance to a neighboring site. It 
has been pointed out by Benioff [Bcnioff 2002 that in a two-dimensional network 
under this restriction, Grover's search [Grover 1996] . [Grover 1997] provides no ad- 
vantage in terms of running time over a classical search due to the intrinsic non- 
locality of Grover's symmetrization. Ambainis et al. proposed a generalized formal- 
ism for quantum walk (QW) based spatial search algorithms and worked out the 
specific case of a two-dimensional cartesian network, obtaining a 0(y/N log N) algo- 
rithm |Kempe et al. 2005] , to which we shall refer to as AKR. Tulsi has proposed an 
improvement to AKR which requires an ancilla qubit and leads to an 0(\/N log N) 
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algorithm in two dimensions [Tulsi 2008 . However, it is not known whether the op- 
timal solution is O(yfN) for the two-dimensional spatial search problem. In contrast, 
it is known that the optimal solution is achieved in higher dimensions, such as 3D- 
grid |Kempe et al. 2005] and the SKW algorithm [Shenvi et al. 200"3] , which searches 
an item within N = 2™ sites arranged in an n-dimensional (n > 2) hypercube. 

There are three ways to cover the plane with regular polygons: squares, hexagons 
and triangles. The resulting regular networks differ in their degree d (number of 
connections per node) which is 4, 3 and 6 respectively. As mentioned before, for 
the rectangular grid (d — 4) the AKR search algorithm finds a marked vertex 
in time 0{sfN log N). A QW- based spatial search algorithm of time complexity 
0(y/N log N) has recently been implemented for a two-dimensional hexagonal net- 
work (d = 3) |Abal et al. 2010] . Both algorithms can be improved to 0(y/ N log N) 
with Tulsi's modification. These results suggest that the degree or connectivity of a 
regular network does not affect the performance of a QW-based search algorithm. 

In this work, a new search algorithm for the case of a triangular network is pro- 
posed and analyzed. The proposed algorithm is a special case of the general framework 
for abstract search |Kempe et al. 2005 [Tulsi 2008] , so we employ the AKR- Tulsi for- 
malism to show that the time complexity of the quantum search on the triangular 
lattice is 0(y/NlogN). This provides further evidence that the degree of the under- 
lying network does not affect the performance of the quantum algorithm. 

The paper is organized as follows. In Section [2] we discuss the implementation 
of a QW on the triangular network. In Section [3] we analyze the time complexity of 
this search algorithm. In Section [4] we perform a numerical analysis of the algorithm. 
Finally in Section [5] we present our conclusions. 



2 QW on a triangular Network 

Let us consider iV sites arranged in a triangular network covering a two-dimensional 
region, as shown in Figure [TJ The network is y/~N x y/~N and periodic boundary 
conditions are assumed. A site on the lattice is located by two integers {m,ri2) 
according to 

r = niai + n 2 a 2 , (I) 

where ai and a 2 are unit vectors forming a 60° angle, as indicated in Figure [TJ These 
integers are such that n.; € [0, y/~N — 1] for i = 1, 2 and thus each of them takes y/N 
different values. 

These sites define an orthonormal set of quantum state vectors, {|ni,ri2)}, which 
span an TV-dimensional Hilbert space Hp. At a given site, there are six possible 
directions of motion which we label with an integer j £ [0, 5], as indicated in Figure [TJ 
The orthonormal states {\j)} span a six-dimensional Hilbert space, He, which we shall 
refer to as the "coin" subspace. The Hilbert space for this problem, % = He ® Hp, 
is 6iV-dimensional. A generic state vector is expressed as 

5 VN-1 

= £ E «S,*\j,*), (2) 

j'=0 m ,n 2 =0 
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Fi gure 1 '. A \/~N X y/~N triangular network (here N = 36) . The N sites form a Bravais lattice and 
the six directions of motion are labelled by an integer j £ [0, 5]. 



where the aj „ are complex amplitudes which satisfy the normalization constraint and 
we have introduced the shorthand notation n = (7x1,7x2). 

The standard QW on this network is implemented with a unitary evolution oper- 
ator of the form 

U = S-(C(g>I P ) (3) 

where S is a shift operator in H (to be specified below), Ip is the identity operation 
in Hp and C is a unitary coin operation in "He- The useful coin operation for spatial 
search problems |Kempe et al. 20051 IShenvi et al. 200"3] is Grover's coin, whose ma- 
trix elements for a ci-dimensional space are = | — 
d = 6, it is given by 



For the particular case 
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Thus we use C — The shift operator implements single-step displacements 

acting on the kets |j; n±, n-i) in the form 



S|0;ni,n 2 ) 
S\l;ni 1 n 2 ) 
S\2;n l7 n 2 ) 
S , |3;ni,n 2 ) 
S , |4;ni,n 2 ) 
5|5;m,n 2 ) 



3;?ii + 1,712)) 
4;m + l,n 2 - 1), 
5;m,n 2 - 1), 
0;ni - 1,712)1 

ljril -1,7X2 + 1), 

2;ni,7x 2 + 1). 



(5) 



Note that S inverts the coin state. This invertion is crucial for the efficiency of the 
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search algorithm described in the next section. Finally, the dynamics of the QW is 
obtained by applying U repeatedly |^(m)) = U m \^f(Q)) for some integer m. 

The standard QW is best analyzed in the Fourier-transformed space. Let us 
consider the reciprocal lattice vectors {gi,g2}, which satisfy the usual requirements 
from condensed matter physics [Kittel 1995) 



A site in this reciprocal lattice is located by 

k = fcigi + fc 2 g 2 , (7) 



with fci, ki integers in [0, v^V — 1]. Let us use the notation k = (fci, k,2) and write a 
generic state vector in the Fourier representation as 

l*> = E4*l**>- (8) 

The kets \ j, k) and \j, n) are related by the discrete Fourier transform 

l*> = ^fE^""!") ( 9 ) 



N ■ 

■xi 

4^J2^ k ) ( 10 ) 



N ■ 

k 

and one can check that (k\h) = uj n ' k /y/~N where u = e 2 '"' 1 /^ and n-k = n\ki + 712^2- 
The action of S on the kets \j, k) of the Fourier representation can be obtained from 
Eqs. © as 

S\Q,k) = w kl \3,k), S\l,k) = w kl - k2 \4,k), 

S\2,k) = w~ k2 \5,k), S\3,k) = w- kl \0,k), (11) 

S\4,k) = w- kl+k2 \l,k), S\5,k) = w k2 \2,k). 

Thus, in the ^-representation S acts diagonally, i.e. S = ^2tSif\k)(k\, where 
Sk = (k\S\k) is the reduction of S to He- Therefore, the evolution operator, Eq. ((3]), is 
also diagonal in the Fourier representation and can be expressed as U = J2% U%\k){k\, 
where Ut — (k\U\k) acts in T-Lc- The matrix elements of the reduced operator can be 
calculated from Eq. fTT) . with the result 

I w - fe i w~ kl w~ kl —2w~ kl w~ kl w~ kl \ 



(12) 
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The characteristic polynomial factors as 

P(A) = (A - 1) 2 (A + 1) 2 (A 2 - 2 cos(0 fc )A + 1), (13) 
where 8k is defined by 

cos(6> fc ) = i ^cos(fci) + cos(fc 2 ) + cos(fc~i - fc 2 )) , (14) 

and hi = ^= for i = 1,2. The eigenvalues are ±1, each with multiplicity 2, and e ±lSk . 
Let us denote by \v±i) and \v±k) the normalized eigenvectors of associated with 
the eigenvalues ±1 and e ±t6k , respectively. Then \v±\, k) and \v±f., k) are eigenvectors 
of U associated with the same eigenvalues, where is defined in Eq. ([9]). 



3 Time complexity of the search algorithm 

We shall use Tulsi's version [Tulsi 2008 of the framework of the abstract search al- 
gorithm |Kempe et al. 2005] to analyze the time complexity of a search algorithm on 
the triangular network. We assume that there is a single marked vertex which 
we want to find. The search algorithm uses a conditional coin operation, which acts 
as —Ic on the searched site |t) and as otherwise. Thus, the modified evolution 
operator is U' = S ■ C, where C acts in H as just described, i.e. 

a = -i c ® i*> <*i + J2 G(6) ® WW- ( 15 ) 

AKR |Kempe et al. 2005] have shown that the evolution of the modified quantum 
walk U' may be analyzed using the eigenspectrum of the standard QW evolution 
operator U. This fact actually reduces the analysis of the search algorithm to a 
tractable eigenproblem for the unitary operator U. 

The evolution operator U' can be written in another form, useful for Tulsi's modi- 
fied algorithm. This modification requires an extra register (an ancilla qubit) used as 
a control for the operators i? t - and U. Using Eq. (|15p . one can show that U' — U ■ i? t - , 
where i? t - = igiv — 2|uc, i)(uc i\, U is given by Eq. ([3]) and \uc) is the uniform su- 
perposition of the computational basis of the coin space. The operators acting on 
the ancilla register are described in Figure [2] where —Z is the negative of Pauli's Z 
operator and 

X S =( C ° S \ Sin " V (16) 
\_ — sin o cos o J y ' 

where cos<5 oc 1/ylogiV. 

The new evolution operator is 

U" = (-Z ® /) • C(U) ■ (Xj ® /) ■ C(Rt) ■ (X s ® /), (17) 

where C(U) and C(Rf) are the controlled operations shown in Figure [2] and / is the 
identity operator in T~L. We will show that U" must be iterated 0(\/N log N) times, 
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Figure 2: Tulsi's circuit diagram for the one-step evolution operator of the quantum 
walk search algorithm. 



taking \ l)\uc)\up) as the initial condition, in order to maximize the overlap with the 
search element. 

The expression for the controlled i? t ~ is C(R^) — I\2N ~ 2j 1, uc,i){l, uc, i\- Let us 
define 



\5i) = Xl\l) = -sin5|0) +cos5|l>, 
then we define a new reflection operator 

R t = (x\ ® /) • C(Rf) ■ (X s ® I) = I 12N -2\S u u c , t)<*i, «c, t\ 
The effective target state is 



(18) 

(19) 

(20) 
(21) 

Note that the eigenspectrum of U is determined from the eigenspectrum of U. In 
fact, for k ^ 



Let us define 



|t> - \Si,u c ,i). 



U= (-Z&I) -C{U). 



U\0)W± k ,k) = -\o)W± k ,k), 
U\l)\v ±k ,k) - e ±w *|l)|i/± fc ,fc). 



(22) 
(23) 



Eigenvectors |Z)|^±i,fc), I = 0, 1 will not be used, because (uc\v±i) = 0. They are 
orthogonal both to the initial condition and to the target. For k — 0, the initial 
condition is the only eigenvector with eigenvalue 1 that will be used. 

The search algorithm consists in applying U" repeatedly taking the initial condi- 
tion as 

|*o> = \l,u c ,u P ) = -= V|l,j,n), (24) 
•'6iV — 



where \uc) and \up) are the uniform superposition in the coin and position spaces 
respectively. The initial condition can be prepared in 0(y/~N) time steps. The number 
of iterations of U" is given by 7r/4a, where a is defined in the following way. The 
eigenvalues of U" that are different from 1 have the form e ±laj , where < ay < n/2. 
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Then a = min{ai , ota, • ■ • }. AKR have shown that it is possible to estimate a knowing 
the eigenspectrum of the evolution operator U. 

The first step in the analysis of the algorithm is to decompose the target state in 
the eigenbasis of U, that is, we have to calculate coefficients ao, at and bk such that 

u c ,t) = -^=|*o} + ^=y"a fe (|l,i/ +fe ,fc} + |l,z/_ fe ,fc}) + 
v jv v jv : — v ' 

fe#0 

^£^|0 > i/+*,*) + |0 > i/_ fc ,fc)). (25) 



. N T 

kjtO 

Coefficients bk are related to the (— l)-eigenspace. Coefficients ak and bk are real. In 
order to satisfy these conditions, eigenvectors \l,v + k,k) for I = 0,1 must be chosen 
appropriately. The procedure is the following. If {uc\v+ k ) has a non-zero phase e lA , 
then \l,v + k, k) must be redefined to e~ lX \l, v+k, k). The same procedure must be 
performed to \l, v-k, k). After those redefinitions, Eq. (|25|) is valid. For the triangular 
network, a straightforward calculation yields 

COS f) 

a k = (l\S 1 )(v ±k \u c ) = ^ (26) 

b k = (0\ Sl ){v ±k \uc) = - S -^ (27) 

This result is remarkably similar to the result obtained for the 2D grid. Note that 
the corresponding expressions for the honeycomb lattice [Abal et al. 2010] are more 
complex, because they depend on k. 

Tulsi has shown that a can be determined using the expression 



(28) 



when N 3> 1. Eqs. (|2"B1 and (|2"T|) show that ak and bk do not depend on k. The 
non-trivial sum inside the above square root may be calculated using 

1 1 N 1 f i-27T-£ ^ 

2 g * is c^yy. dMfci rr^v (29) 

(■27T- 




where we have used -> f // E £ dMfc 2 with e = 2tt^/2/N < 1. Using 

Eq. ([T4| and £ -C 1, Eq. ([29| can be approximated by 



^4 r ^ r~ s . ^ ■■ -w^Ujvbg*. (so) 

16 ^ 7 £ 2 7 e fc 2 + fc 2 - fcifc 2 8 w 

Replacing this result in Eq. (l28l) . using that ao = cos (5, cos 5 = 9 (l/VBgTV) and 
Eqs. ((261) and ([27). we obtain 



- = O (yjVlogiVj . (31) 
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Since the number of iterations of U" is 7r/4a, we conclude that the time complexity 
of the search algorithm is 0(\/N\ogN). It remains to show that the probability to 
find the marked vertex is constant. This is the moment where Tulsi's method shines. 
The overlap of the final state after t = ir/Aa iterations of U" with the target state 

is 

\{5i,u c ,t\(U") t \*o)\=® (minj 1 =,ll| . (32) 

The coefficients bk play no role in this overlap as shown by Tulsi. The calculation 
at this point for the triangular network is similar to one performed by AKR for the 
2D grid, except that ak has the constant factor of cosS. Using that: (1) cot^- < J?, 

k 

(2) i is bounded from below and above by const /(l — cos 6^,) and (3) the result of 

k 

Eq. (1301), we obtain 

| (St , u c , t\ (17")' |* ) | = 6 (min { ^cos 2 (6) log N, l}) . (33) 
Using that cos<5 = 6 (l/v / loglV), we obtain 

\(6 1 ,ucA(U") t \# )\=e(l). (34) 

The probability of finding the marked vertex is constant. We have confirmed this 
result with a numerical simulation in Fig. [4] 



4 Simulation of the search algorithm 

A generic state in the extended space is 

j n 

where a, b are the amplitudes of the state vector extended to include the 2-dimensional 
Hilbert space associated with the ancilla qubit. 

It is straightforward to show that the operator U" in Eq. (fTT]) defines the following 
map: 



*3,n 



a 3fi + ^S A i {a L t sin 2 5 - b i>{ cos 5 sin S) , (36) 



k.l h' 



h,ti' + \^n',i^2 ( a ^,t s i n ^ cos <5 — ^.f cos 2 (S) 



where 5~ | is a Kronecker-delta which selects the searched state. This map allows 
the simulation of the abstract search algorithm in a digital computer. The initial 
condition is taken to be |l)|uc)|up), and the effective target state is \t) = \§i}\uc)\t) , 
as explained in the last section. 
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time step sqrt(NlogN) 

Figure 3: Left panel: probability of finding the searched state against the number of 
time steps for spatial search with (thin line) and without (thick line) Tulsi's modifi- 
cation for N = 400. Right panel: £max against y/N log N using Tulsi's search. Note 
the approximately linear dependence. 



Fortran 90 was chosen as programming language because it provides useful tools 
for large matrix computations and intrinsic functions for complex vector algebra. A 
parallel programming approach using OpenMP was implemented. 

We present the results of the simulations in Figures [3] and |U The left panel in 
Figure |3] shows the time evolution of the probability of finding the searched state, 

i — 1 2 

, both for an algorithm with Tulsi's modification (thin curve), based on U" 
defined in Eq. (|17p. and without it (thick curve), based on U' = S-C with the modified 
coin operator defined in Eq. l[Poj). Note that Tulsi's search yields a smoother curve, 
the first maximum (t max ) of which is reached later but the maximum probability of 
finding the searched element is higher than with the usual spatial search. A change 
in the sign of the local derivative was used to find the maximum point. In the right 
panel of Figure [31 the time at which the maximum probability is reached for Tulsi's 
search is plotted against y/N log N. A straight-line fit to this data has a correlation 
coefficient R 2 = 0.9988. 

In order to test whether the overlap with the searched site in Tulsi's search is con- 

i — 1 2 

stant for large N, in Figure 2] we plot this overlap, \{t\(U") tma,c \l,uc,Up)\ , against 
N. After an initial decay, for N > 2000 the overlap is found to stabilize at approx- 
imately 0.773. These simulations show that Tulsi's search in a triangular network 
works as expected. 

5 Summary and conclusions 

In this work the problem of quantum spatial search on a periodic two-dimensional 
triangular network has been considered for the first time. In this problem, a searched 
item is to be located in a regular triangular network with iV elements. A quantum walk 
operator for this network has been defined and its explicit Fourier representation has 
been found. Its eigenproblem was solved exactly and these results have been used to 
estimate the overlap and running times of the algorithm, according to the generalized 
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Figure 4: Probability of finding the marked vertex | (t\(U / ') tmax |1, uc, up) | against 
N. 



search formalism |Kempe et al. 2005] . This formalism gives a powerful insight on 
the runtime of the algorithm for large values of N. In order to check our analytical 
results and to gain further knowledge on the detailed performance of the spatial 
search, we have implemented its simulation on a classical computer, using standard 
parallel techniques. This allowed us to obtain results for high values of N and see 
how fast the convergence to the theoretical expectations actually is. The simulations 
were implemented both for a modified quantum walk search and for Tulsi's search, 
which uses an extra qubit as a control register. 

Both spatial search algorithms are found to require 0(\/ N log N) steps to reach 
the point at which a measurement yields the searched state with constant probability. 
However, this result is obtained with a constant overlap with the searched state, in 
the case of Tulsi's search. 

Previous work for the spatial search problem in a plane has considered the case 
of a square grid |Kempe et al. 2005] and an hexagonal grid |Abal et al. 2010) . For 
both quantum algorithms the time complexity is 0(y/NlogN). In a sense, this work 
completes the program for the spatial search problem, by providing the details of a 
search in a triangular grid. Since these regular graphs have different degree (d — 4 
for rectangles, d — 3 for hexagons and d = 6 for triangles), this completes the proof 
that the degree of a regular graph does not affect the performance of a spatial search 
algorithm. 

A search algorithm implemented on a real network will have to cope with loss of 
symmetry due to imperfections. The issue of how robust these different spatial search 
protocols are when it comes to searching when a fraction of the links in the network 
is missing is still an important open question which remains for future work. 
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